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INTRODUCTION 


The  Master  Oceanographic  Observation  Data  Set  (MOODS),  maintained  by  the  Naval  Oceano¬ 
graphic  Office  (NAVOCEANO),  Stennis  Space  Center,  Mississippi,  is  the  standard  Navy  database 
for  oceanographic  observations  (references  1  and  2).  Data  contained  in  MOODS  include  observa¬ 
tions  of  thermal  and  chemical  properties  of  the  water  column  from  surface  to  bottom,  the  biota  and 
currents  in  the  water  column,  and  the  underlying  sea  floor  as  collected  by  oceanographic  agencies 
worldwide.  Although  MOODS  provides  the  most  comprehensive  collection  of  oceanographic  data 
available,  these  data  require  analysis  to  provide  insight  into  oceanographic  processes  in  regions  of 
interest  both  to  operational  Naval  forces  and  Navy  Laboratories.  Because  these  regions  frequently 
change  in  response  to  world  politics,  the  Navy  must  be  able  to  examine  regional  ocean  characteristics 
whenever  and  wherever  desired. 

This  document  describes  software  developed  at  the  Naval  Command,  Control  and  Ocean  Surveil¬ 
lance  Center,  RTD&E  Division  (NRaD),  to  process  regional  sets  of  expendable  bathythermograph 
(XBT),  ocean  station  (OSTA),  or  hydrocast  and  conductivity-temperature-depth  (C-T-D)/salinity- 
temperature-depth  (S-T-D)  observations  resident  in  MOODS.  Included  software  offers  a  wide  vari¬ 
ety  of  options  to  filter,  process,  analyze,  and  display  data. 

The  majority  of  programs  described  in  the  document  are  installed  on  a  Sun  SPARC  10  (host  code 
WAHOO);  while  two  routines  are  programmed  for  the  MS/DOS  operating  system  installed  in  a 
386-series  microcomputer.  Operating  guidelines  for  these  routines  are  provided  in  appendix  A  and 
appendix  B,  respectively. 


ARCHIVED  DATA 

The  original  unclassified  MOODS  data  are  located  on  9-track  tapes  located  in  the  NRaD  Code  713 
secure  storage  facility.  Processed  regional  files  are  stored  on  WAHOO  and  may  be  accessed  by 
many  of  the  programs  described  in  this  document.  Criteria  for  the  construction  of  regional  files  were 
based  on  two  factors:  (1)  natural  geographic  divisions  such  as  shorelines,  island  chains,  continental 
shelves,  etc.,  and  (2)  the  number  of  data  available.  Regions  are  not  constrained  to  rectangles,  but 
may  be  downloaded  into  complex  shapes  using  (1)  least-squares  fit  of  latitude/longitude  and  (2)  a 
combination  of  subareas  using  the  MS/DOS  COPY  command.  An  example  of  a  regional  file  is 
SOCAL.MBT,  which  contains  all  available  XBT  data  for  the  Southern  California  Bight.  Figures  1 
through  5  provide  the  names  and  locations  for  MOODS  data  currently  available  on  WAHOO. 

The  MOODS  database  was  revised  in  1993  to  conform  with  security  requirements.  These  data  are 
available  from  NAVOCEANO  Code  N3211  on  diskettes  with  individual  data  files.  These  files  must 
be  preprocessed  on  an  IBM-compatible  personal  computer  using  Routine  NUMDS  (reference  3)  and 
uploaded  to  WAHOO  using  the  Microcopy  Protocol  (MCP)  before  they  can  be  accessed  with  the 
programs  described  in  this  document.^  Theoretically,  larger  data  sets  could  be  transferred  electroni¬ 
cally  from  NAVOCEANO  using  the  File  Transfer  Protocol  (FTP),  but  this  technique  has  not  been 
verified  at  this  time. 

Data  files  saved  by  WAHOO  may  be  either  used  by  other  programs  or  downloaded  via  MCP  to  a 
personal  computer  for  use  with  Program  MDS  (reference  4). 


1  Instructions  for  using  MCP  can  be  obtained  from  the  Custom  Service  Group,  NRaD  Code  02902. 
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DATA  SECURITY 


MOODS  data  contains  both  classified  and  unclassified  data.  However,  because  WAHOO  is  an 
unclassified  system,  only  unclassified  data  can  be  processed  on  this  machine.  Classified  data  may  be 
processed  using  Routines  NUMDS  and  MDS  on  MUNICH,  located  in  the  Bayside  Acoustic  System 
Computing  Center  (BASCC)  in  Building  1,  Bayside.  Permission  to  use  this  system  must  be  certified 
by  the  BASCC  Security  Officer. 


PROGRAM  DESCRIPTION 


ROUTINE  CKREEL 

Routine  CKREEL  examines  a  9-track  reel  and  lists  the  data  included  in  each  10-degree  rectangle 
as  shown  in  figure  6.  This  listing  is  necessary  to  ensure  that  all  data  are  accessed  when  building 
regional  files  using  Routine  XBTO. 

The  program  uses  World  Meteorological  Organization  (WMO)  codes  to  indicate  the  general  posi¬ 
tion  of  the  data.  The  first  digit  of  the  10-degree  code  represents  the  quadrant  of  the  earth  (1:NE, 
3:SE,  5:SW,  7:NW  ),  the  second  the  first  digit  of  latitude,  and  the  last  two  the  first  and  second  digits 
of  longitude.  For  example,  the  first  10-degree  rectangle  shown  in  figure  6  (7215)  is  decoded  as 
20°N,  150°W.  One-degree  WMO  codes  ranging  from  00  to  99  refine  the  position  of  the  first  and  last 
observations  included  in  the  larger  rectangle.  In  this  case,  57  indicates  5  and  7  degrees  for  latitude 
and  longitude,  respectively.  This  adjusts  the  position  of  the  first  observation  to  25  N,  157  W,  which 
is  within  a  few  nautical  miles  of  the  actual  position  of  25°49'N,  157°30'W.  The  listing  also  shows 
that  WMO  rectangle  7215  contains  3081  of  the  10,391  XBT  traces  included  on  reel  U15978. 

ROUTINE  XBTO 

Routine  XBTO  is  the  basic  routine  for  downloading  data  from  9-track  tape  reels  to  regional  files. 
Data  for  the  files  shown  in  figures  1  through  5  were  generated  using  this  program.  The  strength  of 
this  program  is  the  ability  to  use  algorithms  generated  by  Routine  LSTSQRS  to  represent  nonlinear 
boundaries  when  defining  regional  limits.  Figure  7  illustrates  output  from  this  program  listing  the 
number  of  observations  for  each  type  data  downloaded.  File  names  are  the  same  for  all  data  formats, 
with  type  of  observation  indicated  by  file  extension  (.MBT,  .MOS,  or  .MST).  Algorithms  used  to 
define  the  region  are  printed  as  well  as  administrative  information. 

ROUTINE  DWNLD 

Regional  data  files  may  be  edited  for  use  in  the  analysis  routines  described  below  or  for  download¬ 
ing  to  a  personal  computer  for  later  processing.  In  the  latter  case,  smoothing  of  XBT  and  C-T-D/S- 
T-D  data  substantially  reduces  the  time  required  to  download  the  file.  The  limited  numbers  of  data 
points  in  a  hydrocast  (maximum;  32)  negates  the  advantage  of  processing  these  data  with  DWNLD. 

ROUTINE  XBT2A 

Routine  XBT2A  is  the  major  data  analysis  routine  for  XBT  data.  Statistical  output  is  in  two  sec¬ 
tions:  (1)  temperature  at  nine  levels,  sonic  layer  depth^,  and  temperature  gradient  between  200  and 
300  m  (DT23),  and  (2)  characteristics  of  the  near-surface  sound  channel.  Choice  of  periods 
(biweekly  to  annual)  and  geographic  filters  are  available.  Frequency  of  phenomena  such  as  zero 

2  Definitions  for  acoustic  parameters  used  in  the  programs  described  in  the  text  are  given  in  appendix  C. 
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sonic  layer,  half  channel,  and  useful  sound  channels  are  printed  below  the  statistical  information. 
Salinity  from  three  sources  may  be  used  for  sound  speed  computations:  (1)  constant  salinity  of  35 
parts  per  thousand  (o/oo),  (2)  a  single  user-provided  salinity  profile,  or  (3)  seasonal  profiles  gener¬ 
ated  using  Routine  TSZ.  The  distribution  of  observation  per  30-minute  rectangles  over  a  10-  by 
10-degree  area^  may  be  printed  and  saved.  Examples  of  statistical  printout  as  a  function  of  month 
are  shown  as  figures  8  and  9.  Figure  10  shows  a  typical  miniplot  printout.  In  this  plot,  limited  data 
to  the  southwest  (lower  left  corner  of  the  plot)  shows  the  effect  of  limited  observations  owing  to 
shallow  water,  whereas  limited  observations  to  the  northeast  are  the  result  of  a  data  filter.  Crude 
plots  of  temperature  values  at  a  selected  depth  (figure  11)  and  DT23  may  be  printed.  The  data  may 
be  stored  in  two  formats:  (1)  the  statistical  data  for  use  with  routines  XSIG  and  CUMLTV  and  (2) 
number  of  observation  per  30-minute  rectangle  for  use  with  Routine  MRGEXY. 

Note:  Sound  channel  statistics  are  computed  for  that  channel  having  the  lowest  cut¬ 
off  frequency.  In  geographic  regions  where  the  deep  sound  channel  approaches  the 
surface.  Program  MDS  cannot  reliably  discern  the  difference  between  deep  and  sec¬ 
ondary  sound  channels.  This  problem  is  mitigated  by  the  fact  that,  for  most  tactical 
purposes,  sonic  energy  travels  equally  well  in  both  ducts. 

Note:  Because  of  the  considerable  variability  that  occurs  in  the  upper  layer  of  the  sea, 
statistical  computations  should  be  made  for  1 -month  periods  only. 

Note:  Remember  that  values  computed  with  less  than  30  samples  are  not  reliable  sta¬ 
tistically. 

ROUTINE  TSZ 

Hydrocast  and  C-S-T/S-T-D  data  are  used  to  provide  (1)  a  crude  T-S  scatter  plot  at  a  user- 
defined  depth  and  (2)  seasonal  salinity  profiles  with  respect  to  depth.  The  profiles  may  saved  as  a 
source  of  salinity  input  into  Routine  XBT2A.  Examples  of  printout  from  this  program  are  shown  as 
figures  12  and  13. 

ROUTINE  XBT3A 

This  routine  uses  hydrocast  and  C-S-T/S-T-D  data  to  (1)  generate  miniplots  similar  to  those  in 
Routine  XBT2A  and  (2)  mean,  standard  deviation,  and  number  of  observations  per  1 -degree  rectan¬ 
gle  for  temperature,  salinity,  sound  speed,  or  density.  The  miniplot  may  be  stored  for  later  use  by 
Routine  MRGEXY.  An  example  of  output  in  the  latter  configuration  is  shown  in  figure  14. 

ROUTINE  XSIG 

Level  of  significance  between  parameters  from  data  sets  generated  and  saved  by  XBT2A  can  be 
determined  using  Routine  XSIG.  This  can  be  particularly  helpful  in  determining  similarities  and  dif¬ 
ferences  between  data  sets,  thereby  permitting  inferences  about  a  geographic  region.  For  example, 
significance  differences  between  temperature  at  most  levels  from  adjacent  areas  strongly  suggests 
presence  of  an  oceanic  front.  Conversely,  if  no  significance  is  found,  the  data  are  from  the  same 
water  mass  and  may  be  combined  using  Routine  CUMLTV.  This  implies  that  careful  examination  is 
required  within  each  area  to  separate  the  data  into  separate  water  masses.  A  typical  printout  from 
XSIG  is  shown  in  figure  15. 


3.  Henceforth  referred  to  as  a  “miniplot.” 
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ROUTINE  CUMLTV 

Statistical  data  sets  generated  and  saved  using  Routine  XBT2A  can  be  combined  into  a  single  data 
set  and  printed  with  Routine  CUMLTV.  This  function  is  particularly  useful  when  the  number  of 
observations  in  a  given  area  is  insufficient  to  provide  meaningful  statistics  or  where  the  area  of  inter¬ 
est  includes  more  than  one  region.  Output  is  similar  to  that  of  Routine  XBT2A. 

ROUTINE  MRGEXY 

Routine  MRGEXY  can  be  used  to  determine  water  mass  distribution  by  combining  stored  mini¬ 
plots  generated  by  Routines  XBT2A  and  XBT3A.  Up  to  four  water  masses  may  be  included;  how¬ 
ever,  you  must  ensure  that  all  miniplots  for  a  given  water  mass  were  computed  using  identical 
parameters.  Stored  miniplots  are  adjusted  automatically  for  differences  in  geographic  coverage  to 
correspond  with  the  area  designated  by  the  program.  Output  includes  percentage  of  total  number  of 
observations  by  water  mass  by  30-minute  and  1-degree  rectangle.  An  example  showing  the  percent¬ 
age  of  water  typical  of  the  east  central  North  Pacific  Ocean  is  shown  in  figure  16.  Note  that  the  his¬ 
torical  data  included  in  MOODS  span  several  decades;  therefore,  implied  water  mass  boundaries 
represent  only  the  general  position  of  an  ocean  front,  not  the  location  at  a  specific  time.  Further¬ 
more,  data  from  all  seasons  are  used,  and  ephemeral  fronts  may  be  masked  by  this  process.  A  plot 
showing  the  relative  location  of  ocean  features  seaward  of  California  as  generated  using  output  from 
Routine  MRGEXY  is  shown  in  figure  17. 

ROUTINE  LSTSQRS 

LSTSQRS  is  a  PC  routine  for  least-squares  fit  of  data  sets  such  as  latitude  versus  longitude  or  tem¬ 
perature  versus  salinity.  Routines  XBTO,  DWNLD,  and  XBT2A  use  the  algorithms  generated  by 
LSTSQRS  to  determine  the  boundaries  for  a  geographic  filter.  Figures  18  and  19‘^  show  typical  out¬ 
put  for  the  analysis  and  a  copy  of  the  PC  monitor  showing  the  goodness  of  fit,  respectively. 

ROUTINE  XYPLOT 

XYPLOT  is  a  PC  program  that  produces  plots  using  the  output  from  Routine  XBT2A.  The  data 
may  be  used  to  (1)  plot  monthly  mean  and  standard  deviation  of  temperature,  salinity,  and  acoustic 
parameters  versus  depth,  (2)  frequency  of  observations  such  as  zero  sonic  layer  depth  or  useful 
sound  channels,  and  (3)  monthly  mean  of  multiple  parameters  such  as  temperature  at  different 
depths.  Typical  plots  are  shown  as  figures  20  through  22.  These  plots  are  useful  in  illustrating  the 
magnitude  of  temporal  change. 


These  and  other  graphics  captured  from  a  PC  monitor  were  made  using  Deluxe  Paint  II  Enhanced,  a  product  of  Electronic  Arts,  Inc. 
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Figure  1 .  Regional  files  for  the  northeastern  Pacific  Ocean. 
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Figure  2.  Regional  files  for  the  western  Pacific  Ocean. 


6 


LEGEND: 

1  ECHINAl 
5  SCHINAl 
9  SCHINA3 
13  THAI 
17  NBNGL 


2  ECHINA2 
6  SCHINA2 
10  SUNDA 
14  MALAC 
18  SWBNGL 


3  BASHI 
7  EPHILl 
11  SULU 
15  ANDMN 
19  SEBNGL 


4  TONKN 
8  EPHIL2 
12  CELBS 
16  NMOLUC 


Figure  3.  Regional  files  for  Southeast  Asia  and  the  Bay  of  Bengal. 
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8  WARBN2 
11  NEARBN 
14  PAK 
17  GOMAN 


3  ADENW 
6  SSOMALIA 
9  CARBN 
12  NCARBNl 
15  NWARBNl 
18  STRHOR 


Figure  4,  Regional  files  for  the  Arabian  Sea  and  adjacent  areas. 
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LEGEND: 

1  ALBORAN  2  NAFRICAN  3  ALGERIAN  4  SFRANCE 

5  VAL  6  TYRR  7  STRSIC  8  LIBYA 

9  ADRIATIC  10  lONIANl  11  IONIAN2  12  IONIAN3 

13  SIDRA  14  EMED  15  SEMED  16  LEVANT 

17  NEMED  18  AEGNS  19  AEGNN 

Figure  5.  Regional  files  for  the  Mediterranean  Sea. 
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0 

20-  2N 

160-41W 

END: 

99 

29-25N 

169-50W 

4282 

START: 
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PROGRAM  FINISHED 
FORTRAN  STOP 
$ 

Figure  6.  Typical  listing  of  tape  contents  generated  by  Routine  CKREEL. 
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FILE: 

DATE: 

NUMBER  XBT  TRACES  SAVED:  2365 
NUMBER  HYDROCASTS  SAVED:  1039 
NUMBER  STD/CTD  OBS  SAVED:  1043 
NUMBER  WRITE  ERRORS:  0 

NUMBER  REJECTED  OBS:  1621 
START  TIME:  1220.2  END  TIME:  1228.3 


DATA  FROM  REEL 

X1775  ,  FL  2, 

DOWNLOADED  TO 

FILE  CANNW 

ALGORITHM  TERMS 

NORTH 

SOUTH 

EAST 

WEST 

55 . 00000000 

50 . 00000000 

57.74202000 

140.00000000 

0.00000000 

0.00000000 

1.63092500 

0 . 00000000 

0 . 00000000 

0 . 00000000 

-0.04202001 
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0 . 00000000 
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0.00000000 

0 . 00000000 

0.00000306 

0.00000000 

0.00000000 

0 . 00000000 

0.00000020 

0.00000000 

0.00000000 

0 . 00000000 

0.00000000 

0 . 00000000 

NUMBER  OF  OBS  PROCESSED:  19698 

Figure  7.  Inventory  of  regional  data  downloaded  using  Routine  XBTO. 
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Figure  9.  Statistical  computations  in  the  sound  channel  generated  with  Routine  XBT2A. 
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Figure  10.  Example  of  a  miniplot  produced  by  Routine  XBT2A. 
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Figure  12.  A  T-S  scatter  plot  from  Routine  TSZ. 
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Figure  14.  Statistical  printout  of  temperature  by  1 -degree  rectangles  produced  by  Routine  XBT3A. 


SET 

1: 

EASTPAC 

COLD 

[  5  TO 

10  @  150 

M;  SAL: 

EP3COLD] 

SET 

2 : 

EASTPAC 

WARM 

[  10  TO 

22  @  150 

;  SAL: 

EP3WARM  ] 

EPCOLD 

EPWARM 

SURFACE  TEMP 
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vs 
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vs 

§ 

99 

Figure  15.  Printout  showing  level  of  significance  of  a  given  water  mass  produced  by  Routine 
XBT2A. 
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Figure  16.  Frequency  of  occurrence  of  a  given  water  mass  produced  by  Routine  MRGEXY. 
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Figure  17.  Watermass  analysis  derived  from  combined  frequency  plots. 


DATA  SET:  EASTERN  SIDE  OF  ANDAMAN  SEA 
POWER  COEFFICIENT 

0  127.0306 

1  -9.31768 

2  1.178889 

3  -6.733438E-2 

4  1.456577E-3 

X(I)  Y(I)  L.  S.  Y(I)  L.S.  ERROR 


5.00 

102.40 

102.41 

0.008 

7.00 

100.00 

99.97 

-0.026 

9.00 

99.10 

99.13 

0.031 

11.00 

98.90 

98.89 

-0.015 

13.00 

98.80 

98.80 

0.001 

15.00 

99.00 

99.00 

0.001 

RESIDUAL  SUM  OF 

SQUARES : 

0.0000 

Figure  18.  Least-squares  fit  of  a  geographic  boundary  generated  by  Routine  LSTSQRS. 
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APPENDIX  A 

PROGRAMS  AVAILABLE  ON  WAHOO 


WAHOO  is  located  in  the  General  Purpose  Computer  Center  (GPCC),  Room  B229,  Building  1, 
Bayside.  The  system  operators  will  provide  general  assistance,  including  mounting  tape  reels.  Out¬ 
put  will  be  placed  in  bins  alphabetically  by  user  name  immediately  inside  Room  B229.  Questions 
pertaining  to  the  VMS  operating  system  should  be  directed  to  Code  0299. 


Al.  Most  programs  described  in  this  text  are  found  on  the  Sun  SPARC  10  Computer  network 
WAHOO.  Locations  of  specific  programs  described  in  this  document  are: 

Directory  Routines 


XBTO 

XBTO ,  DWNLD,  CHKRL 

XBTIA 

XSIG 

XBT2A 

XBT2A,  XSIG,  CUMLTV 

XBT3A 

XBT3A,  MRGEXY 

TSZ 

TSZ 

Regional  files  are  stored  in  the  following  manner: 


Directory 

Type  Data 

Extension 

No.  Files 

Blocks 

XBTO 

C-T-D/S-T-D  .MST 

98 

125,604 

XBT2A 

XBT 

.MBT 

109 

275,430 

TSZ 

Hydrocast 

.MOS 

99 

72,745 

A2.  WAHOO  is  accessed  using  the  NRaD  TERMINAL^  command  from  a  personal  computer 
tied  into  the  General  (GCB)  system  using  the  following  procedure: 

T22 

NOSC  Terminal  program  vers  1.15  of  1988  Jan  5 
Configured  for  an  IBM— PC  using  COMl : ,  at  9600  baud, 
emulating  a  DEC  VTlOO 

#LO 

#location  8 

#CA 

#call  a50 

#  CALL  COMPLETED  TO  0A50,1 

*  These  commands,  with  instructions,  may  be  obtained  from  the  Customer  Service  Group,  located  in  Building  A33,  Topside. 

2  In  this  document,  computer  input  is  shown  in  boldface. 


A-1 


<return> 

[  security  information  printed  here  ] 

Username:  jdoe^ 

Password:  abcdef 

Last  interactive  login  on  .  .  .  (date,  time) 

TERM  TYPE  (VTIOO)?  <return> 

You  will  be  in  the  root  directory  upon  successfully  logging  into  WAHOO.  The  desired  directory 
must  be  opened  using  standard  UNIX  operating  procedures: 

SET  DEF  [.directory  name] 

Routines  stored  on  WAHOO  are  run  using  the  command: 

RUN  <routine  name> 

Should  you  want  to  run  another  program  not  in  the  current  directory,  you  must  first  return  to  the  root 
directory  before  opening  the  new  directory: 

SET  DEF  [-] 

SET  DEF  [.new  directory] 

A  listing  of  regional  data  files  may  obtained  by  entering  the  proper  director  as  specified  above  and 
entering: 

DIR  *.<file  extension  .MBT,  MOS,  or  MST> 

Files  may  be  copied  from  one  directory  to  another  by  using  the  command: 

COPY  [FISHER. olddir] file. ext  [ FISHER. newdir] 

WAHOO  is  exited  using  the  commands: 

LOGOUT 

ABCDEF  logged  out  at  .  .  . 

<flO> 

A3.  Routine  CKREEL: 

Start  Routine  CKREEL  by  entering  RUN  CKREEL.  You  then  must  enter  the  reel  number  of  the 
9-track  tape  containing  the  raw  MOODS  data.  The  reel  number  applied  by  NAVOCEANO  is  gener¬ 
ally  used,  but  with  a  suffix  indicating  that  it  is  a  copy  of  the  original. 

ENTER  REEL  NUMBER  TO  BE  CHECKED:  12345W 

3  Use  lower  case  when  entering  administrative  information  into  WAHOO. 
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The  reel  will  be  examined  and  the  number  of  each  type  of  observation  and  the  location  given  by 
Marsden  Square  will  be  printed  as  shown  in  figure  1^.  A  data  listing  also  may  be  obtained  by  exer¬ 
cising  the  second  option  below. 

OPTIONS  AVAILABLE: 

(1)  READ  MOODS  DATA 

(2)  LIST  MOODS  DATA 
INDICATE  OPTION  DESIRED:  1 

If  a  listing  is  desired,  the  sequential  position  in  the  file  of  the  first  and  last  observation  desired  must 
be  given.  This  option  is  not  normally  used  because  of  the  extra  time  required  and  excessive  data  out¬ 
put. 

ENTER  #  OF  FIRST  OB  TO  BE  PRINTED  (XXXX):  1 

ENTER  #  OF  LAST  OB  TO  BE  PRINTED  (XXXX):  100 

Errors  encountered  while  reading  the  data  can  be  listed  if  desired;  however,  this  option  normally  is 
not  exercised. 

DO  YOU  WANT  TO  LIST  ERRORS:  N 

The  routine  may  be  recycled  as  often  as  desired.  However,  the  system  operator  must  be  contacted  to 
dismount  the  original  tape  reel  before  routine  CKREEL  can  be  rerun. 

A4.  Routine  XBTO 

Routine  XBTO  is  started  by  entering  RUN  XBTO.  You  must  then  enter  the  name  of  the  9-track  tape 
as  in  the  previous  section. 

ENTER  THE  NAME  OF  THE  ARCHIVED  DATA  TAPE:  12345W 
IS  THE  NAME  OF  THE  DATA  TAPE  CORRECT?  ANS  Y  OR  N : 

The  system  will  then  request  that  the  desired  reel  be  mounted  and  verify  that  this  has  been  done. 
PLEASE  MOUNT  REEL  12345W 
UNIT  10  OPEN  TO  READ  TAPE  12345W 

MOODS  data  files  normally  contain  data  from  three  sources.  You  may  select  to  use  all  or  only  a  part 
of  the  data: 

INDICATE  TYPE  MOODS  DATA  DESIRED: 

(1)  ALL  DATA  (XBT,  OSTA,  AND  STD) 

(2)  XBT  AND  STD  DATA 

4  Figure  numbers  in  this  appendix  refer  to  figures  from  the  main  text. 

5  Verification  statements  requiring  a  Y  or  N  response  are  used  throughout  the  computer  dialogue.  In  the  interest  of  simplicity,  these  statements 
will  not  be  used  henceforth  in  this  text. 
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(3)  XBT  AND  OSTA  DATA 

(4)  OSTA  AND  STD  DATA 

(5)  XBT  DATA  ONLY 

(6)  OSTA  DATA  ONLY 

(7)  STD  DATA  ONLY 
ENTER  CHOICE  1 

The  data  files  normally  are  contained  in  the  first  file  on  the  tape  reel.  Because  exceptions  do  occur, 
you  will  be  asked  to  specify  the  file  containing  the  desired  data. 

IS  DESIRED  DATA  IN  1ST  FILE  OF  TAPE?  N 

ENTER  FILE  NUMBER  DESIRED:  2 

It  is  expected  that  you  will  want  to  break  the  data  into  numerous  subfiles.  Thus,  you  are  asked  to 
provide  a  file  name  for  each  downloaded  file.  It  is  recommended  that  descriptive  names,  such  as 
geographic  names,  be  used. 

ENTER  NAME  OF  OUTPUT  FILE:  NE_ALSKA 

Separation  of  large  files  into  regional  files  is  the  most  effective  means  of  generating  files  of  man^ 
ageable  size.  This  requires  that  file  boundaries  be  established,  including  geographic  quadrant  of 
the  earth  and  geographic  limits  for  each  of  four  sides. 

Note:  A  region  can  include  data  from  a  single  geographic  quadrant  only.  Should 
a  region  require  data  spanning  the  Greenwich  Meridian,  International  Date  Line, 
or  Equator,  separate  files  should  be  generated  and  merged  later  using  the  MS/ 

DOS  COPY  command.  This  process  also  is  used  to  combine  data  from  different 
tape  reels  or  files. 


INDICATE  QUADRANT: 


(1) 

NORTH 

LATITUDE , 

EAST 

LONGITUDE 

(NE) 

(2) 

SOUTH 

LATITUDE, 

EAST 

LONGITUDE 

(SE) 

(3) 

SOUTH 

LATITUDE , 

WEST 

LONGITUDE 

(SW) 

(4) 

NORTH 

LATITUDE , 

WEST 

LONGITUDE 

(NW) 

ENTER 

NUMBER  OF  DESIRED  QUADRANT:  4 

You  will  be  asked  to  define  each  of  four  boundaries  (upper,  lower,  right,  left)  for  the  region,  with  up 
to  eight  terms  for  each  polynomial.^  For  example: 

BE  SURE  TO  INSERT  DECIMAL  POINT  IN  TERM. 

6  See  Routine  LSTSORS  in  appendix  B  to  generate  polynomial  equations. 
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ENTER  NUMBER  OF  TERMS  FOR  UPPER  BOUNDARY:  1 
ENTER  TERM  1:  55. 

A  grid  with  the  terms  for  each  of  the  four  sides  will  be  shown  for  verification.  If  a  term  is  in  error,  it 
may  be  replaced  by  entering  the  column  (horizontal),  row  (vertical)  position  of  the  erroneous  data 
along  with  the  correct  polynomial. 

DO  YOU  WANT  TO  DOWNLOAD  ADDITIONAL  DATA?  N 

You  can  either  rerun  additional  programs  or  return  to  the  PC  by  using  the  procedure  outlined  in  sec¬ 
tion  A2. 

In  some  instances,  all  data  for  a  desired  region  are  not  included  on  the  same  9-track  tape.  If  this  is 
this  case,  rerun  Routine  XBTO  with  all  parameters  unchanged  except  for  the  names  of  the  regional 
tape  reel  and  output  files.  It  is  suggested  that  the  same  regional  name  be  used,  but  with  numerical  or 
letter  suffixes.  The  file  can  then  be  combined  using  the  UNIX  concatenation  command: 

COPY  /  CONCATENATE  areal . ext , area2 . ext , areaS . ext  area. ext 

Note:  Be  sure  to  specify  the  extension  to  ensure  that  only  the  same  type  of  data 
are  included  in  the  same  file. 


A5.  Routine  DWNLD 

Routine  DWNLD  can  be  run  only  after  regional  files  have  been  constructed  using  Routine  XBTO. 
You  may  modify  file  size,  time  period,  and  water  stratum  as  desired.  It  is  expected  that  the  files  gen¬ 
erated  will  be  downloaded  to  a  personal  computer  using  the  MCP  command;  therefore,  it  is  recom¬ 
mended  that  the  files  be  kept  small  both  to  permit  rapid  downloading  and  easily  managed  files. 

Begin  the  program  by  entering  RUN  DWNLD.  You  must  then  indicate  the  type  data  to  be  processed 
for  downloading: 

SELECT  TYPE  RUN: 

(1)  XBT  TRACES 

(2)  S-T-D  CASTS 

(3)  NANSEN  CASTS 
ENTER  CHOICE:  1 

The  number  of  observations  downloaded  depends  on  the  complexity  of  the  region  and  the  num¬ 
ber  of  observations  available.  In  a  highly  variable  area,  it  is  suggested  that  100  to  200  hundred 
observations  per  month  should  be  adequate  to  represent  normal  conditions. 

NUMBER  OF  OBS  SET  TO  500.  IS  THAT  OK?  N 
ENTER  NUMBER  OF  OBS  DESIRED:  2000 
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Although  the  geographic  area  cannot  be  expanded  from  the  parent  file,  a  data  filter  can  be  used  to 
diminish  the  size  of  the  region  downloaded.  If  this  option  is  selected,  you  will  be  asked  to  enter  four 
geographic  boundaries  as  discussed  in  section  A4. 

DO  YOU  WANT  TO  USE  A  REGIONAL  DATA  FILTER?  Y 

ENTER  NAME  OF  STORED  DATA  FILE;  SWSOJ 

ENTER  NAME  OF  NEW  DATA  FILE  ( 8-LTR  MAX):  SWSOJl 

DATA  WILL  BE  SAVED  IN  FILE  SWSOJl.  IS  THAT  OK?  Y 

ENTER  MINIMUM  DEPTH  DESIRED:  50 

ENTER  MAXIMUM  DEPTH  DESIRED:  800 

DO  YOU  WANT  TO  PROCESS: 


(1) 

ALL  DATA 

(2) 

BY  MONTH 

(3) 

BY  SEASON 

ENTER 

CHOICE: 

3 

ENTER 

DESIRED 

SEASON : 

(1) 

WINTER 

( JAN-MAR) 

(3) 

SUMMER 

( JUL-SEP ) 

(2) 

SPRING 

(APR-JUN) 

(4) 

AUTUMN 

( OCT-DEC ) 

1 

DO  YOU  WANT  TO  RUN  THE  PROGRAM  AGAIN?  N 

You  can  either  rerun  additional  programs  or  return  to  the  PC  by  using  the  procedure  outlined  in  sec¬ 
tion  A2. 

A6.  Routine  XBT2A 

Start  Routine  XBT2A  by  entering  RUN  XBT2A.  You  must  then  enter  a  descriptive  title  for  the 
output.  Each  title  should  describe  the  region  and  physical  characteristics'^  of  the  run.  Specifics,  such 
as  name  and  type  data  contained  in  the  data  file,  watermass  characteristics,  and  geographic  filters 
selected,  etc.,  are  included  on  the  last  sheet  of  the  printer  output.  However,  in  the  likelihood  that  sev¬ 
eral  runs  will  be  used  for  the  same  data  set,  it  is  suggested  that  relevant  information  be  included  in 
the  title.  A  maximum  of  40  characters  are  permitted  in  the  title. 

ENTER  TITLE  OF  OUTPUT;  SO  CAL  DEEP,  T200;  5  TO  12;  HIST  SAL 

Filters  are  required  for  each  run.  You  may  use  broad  filters  to  include  all  data  or  may  limit  the  data 
examined  by  entering  specific  limits. 

^  Acoustic  parameters  generated  by  this  program  are  defined  in  appendix  B. 
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The  optimal  period  for  statistical  processing  is  on  a  monthly  basis.  Longer  periods  should  be  used 
only  to  determine  the  general  parameter  minimas  and  maximas.  The  biweekly  option,  which  uses  24 
15.2-day  periods,  can  be  used  to  determine  trends  over  short  periods  in  areas  where  observations  are 
plentiful. 


ENTER 

STATISTICAL  PERIOD  DESIRED: 

(1) 

ANNUAL 

(2) 

SEMIANNUAL 

(3) 

QUARTERLY 

(4) 

BIMONTHLY 

(5) 

MONTHLY 

(6) 

BIWEEKLY  5 

USE  METRIC  UNITS  FOR  ALL  FUTURE  INPUT 

You  may  select  the  maximum  depth  (MAXD)  to  be  examined  in  the  statistical  output.  This  selection 
is  not  trivial,  because  MAXD  also  determines  the  minimum  depth  required  for  included  observa¬ 
tions,  For  example,  an  observation  extending  to  356  m  would  be  accepted  if  MAXD  were  300  m, 
but  not  if  the  layer  were  400  m.  This  feature  can  be  used  to  assist  separation  of  deep-water  from 
shallow-water  observations. 

ENTER  DEEPEST  STATISTICAL  LAYER  TO  BE  EXAMINED 

SHOULD  BE  50,  100,  150,  200,  300,  400  M:  400 

It  may  be  desirable  to  separate  data  as  a  function  of  water  mass.  This  is  done  by  selecting  a  depth 
and  temperature  limits  for  the  filter.  For  example,  a  temperature  range  of -2°  to  30°  C  at  a  filter 
depth  of  200  m  would  ensure  that  all  XBT  traces  from  the  area  would  be  extracted  in  temperate 
regions  of  the  ocean.  Filter  depth  should  be  deep  enough  to  be  free  of  significant  annual  temperature 
changes,  but  shallow  enough  to  accept  most  of  the  data.  Normally  a  depth  of  150  or  200  m  fulfills 
these  criteria.  Should  only  observations  collected  from  the  continental  slope  and  adjacent  shelf  be 
desired,  you  might  specify  a  minimum  and  maximum  depths  of  50  m  and  300  m,  respectively.  In 
this  case,  truncated  observations  collected  in  deep  water  could  be  eliminated  by  selection  of  a  geo¬ 
graphic  filter.  Values  falling  within  this  filter  will  be  accepted,  all  others  discarded.  If  water  mass 
values  are  not  available  from  the  physical  oceanographic  literature,  the  scatter  diagram  option  of 
Routine  TSZ  in  section  A7  is  used  to  approximate  the  temperature  limits. 

ENTER  DEPTH  OF  TEMPERATURE  FILTER;  200 

ENTER  MIN  TEMP  FOR  FILTER:  5 

ENTER  MAX  TEMP  FOR  FILTER:  12 

The  program  examines  the  temperature  gradient  between  200  and  300  m  to  permit  you  to  separate 
different  water  masses  with  similar  temperature  values  at  the  same  teat  depth.  This  feature  is  useful 
when  neighboring  water  masses  have  similar  criteria  at  filter  depth,  but  not  at  deeper  levels. 
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DO  YOU  WANT  TO  CHANGE  THE  MIN  AND  MAX  TEMP  GRADIENTS 

FROM  -10.0  TO  1.0  ?  Y 

ENTER  MIN  TEMP  GRAD  FOR  FILTER:  -6 

ENTER  MAX  TEMP  GRAD  FOR  FILTER:  0 

The  name  of  the  data  file  downloaded  using  Routine  XBTO  must  be  entered.  The  file  extension 
.MBT  is  not  required.  If  the  file  cannot  be  found,  it  either  does  not  exist  or  has  not  been  copied  from 
directory  XBTO  to  directory  XBT2A. 

ENTER  NAME  OF  DATA  FILE;  SOCAL 

The  maximum  depth  criteria  (ZMAX)  permits  you  to  eliminate  data  that  extend  below  a  given  depth. 
For  example,  the  default  value  of  800  m  would  accept  data  from  T-4  (maximum  depth;  460  m)  and 
T-7  (760  m)  probes,  but  would  not  include  data  from  T-6  (2500  m)  probes.  This  option  can  be  use 
to  separate  shallow  water  data  from  deep  water  data.  However,  it  can  also  be  used  m  deep  water  to 
maximize  the  number  of  probes  available,  for  example,  of  a  near-surface  feature  such  as  layer  depth 
in  Northern  Hemisphere  summer.  ZMAX  must  be  equal  or  greater  than  MAXD;  if  not,  it  will  be 
recalculated  to  MAXD  plus  100  m  unless  directed  otherwise. 

MAX  DEPTH  IS  SET  AT  800  M.  IS  THIS  OK?  Y 

The  input  is  then  displayed  and  you  are  given  the  opportunity  to  make  changes.  Proper  use  of  this 
option  can  save  considerable  run  time.  When  Routine  XBT2A  is  recycled,  this  menu  will  appear 
with  the  parameters  from  the  previous  run,  thus  saving  the  effort  of  re-entering  each  number.  For 
example,  if  a  second  water  mass  is  to  be  examined,  only  the  minimum  and  maximum  values  of  the 
temperature  filter  and  possibly  the  values  for  the  temperature  gradient  filter  need  modification. 

PERIOD  =  1  MOS. 

OBS  LESS  THAN  400  M  NOT  ACCEPTED 
TEMP  FILTER  DEPTH  =  200  M 

NO  OBS  DEEPER  THAN  800  M  WILL  BE  EXAMINED 
MIN  AND  MAX  TEMP  AT  TEST  DEPTH  ARE  5.00  AND  12.00 
MIN  AND  MAX  TEMP  GRADIENTS  BTWN  200  AND  300  M  ARE  -6.0 
AND  0.0 

DATA  FILE  IS  SOCAL 

ARE  THE  INPUT  DATA  CORRECT?  ANS  Y  OR  N  Y 

The  statistical  data  for  each  period  may  be  saved  for  later  processing  by  Routines  CUMLTV  or 
XSIG.  If  the  data  are  stored,  you  will  be  asked  to  enter  a  name  for  the  statistical  file.  The  file  exten¬ 
sion  .STA  will  appended  to  the  file  name  automatically. 
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DO  YOU  WANT  TO  SAVE  THE  STAT  DATA?  Y 

ENTER  NAME  OF  STAT  FILE:  SCALDP 

A  10-degree  by  10-degree  miniplot  showing  the  number  of  observations  per  half  degree  may  be 
printed  if  desired.  Note  that  the  data  cannot  be  used  where  quadrants  of  the  earth  intersect.  The  pro¬ 
gram  will  generate  the  proper  geographic  coordinates  for  the  upper  right  corner  of  the  plot  upon 
entry  of  data  defining  the  lower  left  corner. 

These  data  may  be  saved  for  use  in  Routine  XBT3A.  If  the  grid  is  saved,  it  is  recommended  that  the 
same  name  as  used  for  stored  statistical  data  be  used  to  keep  you  from  being  overwhelmed  by  differ¬ 
ent  file  names.  The  extension  .XPL  will  be  added  automatically  for  a  stored  XBT  miniplot. 

DO  YOU  WANT  A  10X10  MINIPLOT?  Y 

MINIPLOT  WILL  NOT  WORK  WHERE  EQUATOR,  GREENWICH  MERIDIAN, 

OR  INTERNATIONAL  DATE  LINE  INTERSECTS  PLOT. 

ENTER  LAT  OF  LWR  LEFT  CORNER  IN  DEGREES  (EG,  35):  30 

ENTER  HEMISPHERE  OF  LL  CORNER  (  N  OR  S  ) :  N 

ENTER  LON  OF  LWR  LEFT  CORNER  IN  DEGREES  (EG,  125):  40 

ENTER  HEMISPHERE  OF  LL  CORNER  (  E  OR  W  ) :  W 

DO  YOU  WANT  TO  SAVE  THE  MINIPLOT  TO  THE  DISK?  Y 

ENTER  NAME  OF  MINIPLOT  FILE  (10-LTR  MAX):  SCALDP 

Histograms  of  the  distribution  of  temperature  at  test  depth  and  the  temperature  gradient  between  200- 
and  300-m  levels  may  be  generated.  Both  of  these  plots  are  useful  in  determining  the  presence  of 
more  than  one  water  mass.  Indications  of  the  presence  of  multiple  water  masses  are  (1)  multiple 
maxima  and  minima  in  lieu  of  a  bell-shaped  curve,  and  (2)  a  wide  range  of  values  (e.g.,  a  10-degree 
temperature  spread  without  a  noticeable  maximum. 

DO  YOU  WANT  TO  PRINT  A  TEMPERATURE  HISTOGRAM  AT  DEPTH?  Y 

DO  YOU  WANT  TO  PRINT  A  TEMP  GRADIENT  HISTOGRAM  AT  DEPTH?  Y 

All  acoustically  related  parameters  can  be  generated  using  historical  salinity  values  in  lieu  of 
default  salinity  of  35  parts  per  thousand  (o/oo).  Historical  profiles  may  be  entered  by  either  of 
two  methods:  (1)  quarterly  profiles  at  standard  oceanographic  depths  as  imported  from  Routine 
TSZ  or  (2)  a  single  profile  entered  manually.  Salinity  will  be  interpreted  for  each  depth  incre¬ 
ment  in  the  XBT  profile  regardless  of  the  method  used. 

DO  YOU  WANT  TO  USE  A  HISTORICAL  SALINITY  PROFILE?  Y 

SELECT  TYPE  INPUT  DESIRED: 

( 1 )  RECALL  STORED  PROFILE 

( 2 )  KEYBOARD  ENTRY  OF  PROFILE 
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ENTER  CHOICE:  1 

The  most  common  means  of  salinity  input  is  using  ocean  station  data  downloaded  by  Routine  TSZ. 

If  this  option  is  selected,  then  the  file  (with  extension  .SAL)  first  must  be  copied  from  Directory  TSZ 
using  the  COPY  command  described  in  section  A2. 

DO  YOU  WANT  TO  USE  A  REGIONAL  DATA  FILTER?  Y 

ENTER  NAME  OF  SALINITY  FILE:  SOCAL 

If  a  manually  entered  salinity  profile  is  desired,  you  will  be  requested  to  enter  salinity  values  at  the 
sea  surface,  50, 100,  150,  200,  300,  400,  500,  600,  and  800  m.  For  example: 

ENTER  SALINITY  AT  0  M:  33.256 

ENTER  SALINITY  AT  50  M:  - 

The  program  provides  the  option  of  excluding  data  outside  a  set  of  four  geographic  boundaries.  This 
procedure  is  useful  in  limiting  the  data  set  to  a  particular  area  such  as  the  continental  shelf  or  on  one 
side  of  the  historical  location  of  an  oceanic  front  to  ensure  that  most  samples  are  representative  of  a 
particular  water  mass.  The  procedure  used  to  apply  a  geographic  filter  in  Routine  XBTO  is  defined 
in  section  A4. 

DO  YOU  WANT  TO  USE  A  REGIONAL  DATA  FILTER:  N 

In  some  instances,  you  might  want  to  download  temperature  data  from  a  specific  depth  for  further 
analysis.  For  example,  data  downloaded  from  this  program  would  be  useful  to  examine  surface  tem¬ 
perature  using  harmonic  analysis. 

DO  YOU  WANT  TO  DOWNLOAD  TEMPERATURE  AT  A  SPECIFIC  DEPTH?  Y 
ENTER  NAME  FOR  TEMPERATURE  FILE:  SCALDP2 
ENTER  DEPTH  FOR  TEMP  DATA:  200 

The  statistical  listings  normally  are  printed  on  the  system  standard  (SYSSPRINT)  using  the  day¬ 
time-group  at  the  time  designated,  plus  a  system  identifier  (example:  073030$PRINT.LIS).  Up  to 
five  copies  of  the  listings  may  be  printed;  the  default  value  is  a  single  copy. 

DO  YOU  WANT  OUTPUT  SENT  TO  A  PRINTER?  Y 

ENTER  NUMBER  OF  COPIES  (MAX  =5):  1 

You  have  the  option  to  make  as  many  runs  as  desired.  If  additional  runs  are  made,  an  abbreviated 
data  menu  is  used  to  set  the  program  filters.  This  provides  a  considerable  savings  in  run  set  up  time. 

DO  YOU  WANT  TO  RUN  THE  PROGRAM  AGAIN?  Y 
A7.  Routine  TSZ 

The  initial  set  up  procedure  for  Routine  TSZ  is  similar  to  that  for  Routine  XBT2A  (section  A6 
above)  except  that  only  hydrocast  (Nansen  cast,  file  extension  .MOS)  and  C-T-D/S-T-D  (.MST) 
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data  can  be  used  as  input.  The  ability  to  select  maximum  and  minimum  parameter  values,  a  regional 
filter,  period,  etc.,  permits  you  to  customize  the  run  to  provide  data  for  a  particular  time  and  place. 

After  entering  the  plot  title,  you  will  be  asked  to  indicate  the  type  data  to  be  processed: 

SELECT  TYPE  DATA  DESIRED: 

(1)  HYDROCAST 

(2)  C-T-D/S-T-D 
ENTER  CHOICE:  1 

You  must  then  designate  the  type  output  desired. 

INDICATE  TYPE  PLOT  DESIRED: 

(1)  TS  SCATTER  PLOT 

(2 )  TS  PROFILE (S) :  1 

The  scatter  plot  provides  a  rough  estimate  of  the  relation  between  salinity  (horizontal  plot  axis)  and 
temperature  (vertical  axis)  at  a  specified  depth  as  shown  in  figure  12.  Note  that  the  plot  is  limited  to 
a  20  X  20  array,  with  the  printed  number  representing  the  number  of  observations  falling  within  each 
bin.  Thus,  the  granularity  of  each  bin  is  a  function  of  the  temperature  and  salinity  ranges  selected. 
For  example,  minimum  and  maximum  temperature  values  set  at  0°  and  20°C  would  result  in  each 
temperature  bin  having  a  value  of  1  degree;  whereas  settings  of  0°  to  10°C  would  have  0.5-degree 
bins.  Thus,  when  determining  water  masses  delineators,  it  might  be  necessary  to  initially  plot  the 
entire  spectrum  of  TS  ranges  followed  by  one  or  more  plots  using  finer  scales  such  as  0.25-degree 
bins.  If  a  finer  plot  is  desired,  you  should  download  S-T-D  data  to  a  PC  and  run  Routine  MDS  as 
described  in  reference  4. 

If  the  TS  Profiles  option  is  selected,  the  number  of  observations,  mean  and  standard  deviation  are 
plotted  at  each  standard  depth  for  temperature,  salinity,  and  sound  speed  for  each  period.  The  pro¬ 
files  also  provide  a  general  indication  of  the  axial  depth  of  the  deep  sound  channel.  Salinity  profiles 
generated  with  this  option  may  be  saved  for  use  with  Routine  XBT2A  by  answering  affirmatively  at 
the  prompt: 

DO  YOU  WANT  TO  SAVE  THE  SALINITY  PROFILES:  Y 
ENTER  THE  NAME  OF  THE  SALINITY  FILE:  SOCAL 
Program  termination  is  similar  to  that  used  for  Routine  XBT2A  in  the  previous  section. 


A8.  Routine  XBT3A 

Routine  XBT3A  generates  a  10-  by  10-degree  plot  similar  to  that  produced  by  Routine  XBT2A 
with  two  additions:  (1)  mean  values  and  standard  deviation  of  temperature,  sound  speed,  salinity,  and 
sigma-t  plus  number  of  observation  are  plotted  as  a  function  of  1 -degree  rectangles  and  (2)  all  data 
files  (.MBT,  .MST,  and  .MOS)  may  be  accessed.  Miniplots  may  be  printed  and  saved  as  in  Routine 
XBT2A  (section  A6)  for  later  use  in  Routine  MRGEXY  (Section  All).  Miniplots  saved  with  this 
program  have  the  extension  .XPL  of  XBT  data,  .OPL  for  hydrocast  data,  and  .SPL  for  C-T-D/S- 
T-D  data. 
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The  program  is  initiated  by  entering  RUN  XBT3A.  The  set  up  for  Routine  XBT3A  is  similar  to 
Routine  XBT2A,  with  the  exception  that  the  type  data  and  output  must  be  identified. 

SELECT  TYPE  DATA  DESIRED: 

(1)  XBT 

( 2 )  HYDROCAST 

(3)  S-T-D/C-T-D 
ENTER  CHOICE:  2 

If  XBT  data  are  selected,  you  have  the  option  of  printing  either  statistical  data  or  number  of  observa¬ 
tions  at  one  of  eight  depths  (0,  25,  50, 100, 150,  250,  300,  or  400  m).  If  hydrocast  or  S-T-D/C-T-D 
data  are  selected,  you  may  process  additional  information: 

SELECT  TYPE  DATA  TO  BE  PROCESSED: 

(1)  TEMPERATURE 

(2)  SOUND  SPEED 

(3)  SALINITY 

(4)  SIGMA-T 

(5)  #  OBS  PER  30-MINUTE  RECTANGLE 
ENTER  CHOICE;  2 

A9.  Routine  XSIG 

Routine  XSIG  uses  Student’s  t-test  to  determine  significance  between  data  sets  generated  with 
Routine  XBT2A  and  saved  in  .STA  files.  Start  the  program  by  entering  RUN  XSIG.  You  then  will 
be  asked  to  enter  the  names  of  the  two  files  that  will  be  used  in  the  comparison,  e.g., 

ENTER  NAME  OF  FIRST  STORED  STAT  FILE;  EPAC3 

ENTER  NAME  OF  SECOND  STORED  STAT  FILE;  EPAC3A 

A  list  of  temperature  and  acoustic  parameters  available  for  comparison  will  appear  on  the  screen  and 
the  desired  choice  indicated. 

PARAMETERS  AVAILABLE  FOR  PROCESSING; 


(1) 

SEA  SURFACE 

TEMPERATURE 

(2) 

TEMPERATURE 

AT 

25 

M 

(3) 

TEMPERATURE 

AT 

50 

M 

A-12 


(  )  SONIC  LAYER  DEPTH 
(  )  TEMP  AT  SONIC  LAYER  DEPTH 
(  )  TEMP  GRAD  BTWN  200  AND  300  M 
(  )  CUT-OFF  FREQUENCY  IN  SURFACE  LAYER 
(  )  SECONDARY  SOUND  CHANNEL  AXIS  DEPTH 
(  )  TEMPERATURE  AT  AXIAL  DEPTH 
(  )  SOUND  CHANNEL  THICKNESS 
(  )  SOUND  CHANNEL  MAGNITUDE 
(  )  SOUND  CHANNEL  CUT-OFF  FREQUENCY 
ENTER  SELECTION:  2 

The  program  may  be  reiterated  as  often  as  necessary  without  renaming  the  data  files.  However,  you 
have  the  option  of  changing  the  data  files  and  making  addition  runs  if  desired. 

AlO.  Routine  CUMLTV 

Routine  CUMLTV  combines  statistical  data  saved  previously  with  Routine  XBT2A.  The  com¬ 
bined  data  is  saved  under  a  new  file  name  with  the  extension  .STA.  The  maximum  number  of  data 
sets  that  can  be  combined  with  this  routine  is  nine. 

After  starting  the  program  by  entering  RUN  CUMLTV,  the  following  information  is  required: 

ENTER  TITLE  OF  COMBINED  STATISTICAL  SETS: 

EPAC  COMPOSITE 

ENTER  NAME  OF  NEW  FILE:  EAST  PAC 

ENTER  NAME  OF  FIRST  STORED  STAT  FILE  1:  EPACl 


The  files  will  be  combined  and  a  statistical  output  will  be  printed  for  the  combined  files  using  a  tabu¬ 
lar  output  similar  to  that  used  for  Routine  XBT2A. 

All.  Routine  MRGEXY 

Water  mass  distribution  can  be  determined  using  Routine  MRGEXY  by  combining  miniplots  gen¬ 
erated  by  Routines  XBT2A  and  XBT3A  into  a  single  set  of  plots.  After  starting  the  program  by 
entering  RUN  MRGEXY,  the  set  up  is  similar  to  the  other  programs: 
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ENTER  TITLE:  CALIFORNIA  CURRENT 

ENTER  POSITIONAL  INFO  OF  BASE  PLOT: 

ENTER  LAT  OF  LWR  LEFT  PLOT  CORNER:  30 
ENTER  HEMISPHERE  OF  LAT:  N 
ENTER  LON  OF  LWR  LEFT  PLOT  CORNER:  125 
ENTER  HEMISPHERE  OF  LON:  W 
ARE  LAT  AND  LON  CORRECT?  Y 

Note  that  the  geographic  location  of  the  composite  chart  need  not  be  the  same  as  that  of  stored  mini¬ 
plots  called  by  this  program;  an  offset  algorithm  adjusts  stored  data  to  fit  the  desired  coordinates. 

All  data  for  a  single  water  mass  must  be  input  before  the  next  water  mass  can  be  addressed.  After 
entering  the  water  mass  title,  you  must  enter  the  name  and  designate  the  source  (.XPL,  .OPL,  or 
.SPL)  of  each  stored  miniplot  file  sequentially.  A  data  set  created  from  XBT  data  (file  extension: 
.XPL)  is  in  the  same  format  as  a  set  generated  from  hydrocast  (.OPL)  or  C-T-D/S-T-D  (.SPL)  data. 
Thus,  these  sets  may  be  summed  to  determine  the  total  number  of  observations  passing  the  desig¬ 
nated  filter.  However,  you  must  first  copy  the  data  sets  into  Directory  XBT3A  using  the  COPY  com¬ 
mand  (section  A2). 

ENTER  ID  FOR  WATER  MASS  #1:  EPAC 

ENTER  NAME  OF  STORED  DATA  FILE  1  FROM  EPAC:  EPAC3 
SELECT  TYPE  DATA  DESIRED: 


(1) 

XBT 

(2) 

OSTA 

(3) 

STD 

ENTER 

CHOICE  1 

After  each  data  set  has  been  read,  you  will  be  asked  if  there  are  additional  files  for  this  water  mass. 
Answer  affirmatively  if  additional  files  are  to  be  included;  negatively  if  all  files  have  been  read. 

ARE  THERE  MORE  FILES  FOR  THIS  WATERMASS?  N 

If  all  data  have  been  read  for  a  given  water  mass,  you  must  indicate  whether  data  from  another  water 
mass  is  to  be  processed. 

ARE  THERE  ADDITIONAL  WATER  MASSES?  N 

The  data  will  then  be  combined  and  miniplots  printed  showing  the  total  number  of  observations  and 
the  percentage  of  each  water  mass  by  30-minute  and  1 -degree  rectangles. 
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APPENDIX  B 

SUPPLEMENTAL  PC  ROUTINES 


The  routines  discussed  in  this  appendix  are  written  in  PowerBASIC^  for  use  on  an  IBM-compat¬ 
ible  personal  computer.  Output  is  formatted  for  a  LaserJet  Series  II  printer. 


Bl.  Routine  LSTSORS 

Routine  LSTSQRS  is  activated  by  entering  LSTSQRS.  You  will  be  asked  to  enter  the  plot  title 
followed  by  the  input  data. 

ENTER  NAME  OF  DATA  SET:  E.  ANDAMAN  SEA 

Be  sure  that  the  data  are  entered  in  the  proper  format;  if  you  want  latitude  as  a  function  of  longitude, 
you  must  first  enter  longitude,  then  latitude.  After  the  data  have  been  accepted,  the  data  will  be 
sorted  in  ascending  order  based  on  the  magnitude  of  the  x  coordinate. 

REGRESSION  IS  Y=F(X)  FORMAT. 

ENTER  X  AND  Y  SEPARATED  BY  A  COMMA;  50  PTS  MAX.  USE  -999,0 
TO  END. 

1  5,102.4 

2  7,100.1 


n  -999,0 

After  all  the  data  have  been  entered,  a  grid  of  the  inputs  will  be  printed  on  the  screen.  You  can  either 
accept  or  modify  the  data  at  that  time. 

AVAILABLE  OPTIONS: 

(1)  ACCEPT  DATA 

(2)  CORRECT  DATA 

(3)  ADD  A  POINT 

(4)  DELETE  A  POINT 

ENTER  NUMBER  OF  DESIRED  OPTION:  1 


'  A  product  of  Borland  International. 
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If  a  point  requires  correction,  enter  the  integer  of  the  point  in  error  and  correct  values: 

ENTER  INTEGER  OF  INCORRECT  ENTRY  AND  CORRECTED  DATA: 

3,9,99.1 

If  a  new  data  point  is  desired,  enter  the  new  coordinates  and  the  program  will  automatically  place 
them  in  the  proper  order  in  the  array  based  on  the  magnitude  of  the  term  x. 

ENTER  NEW  X  AND  Y  VALUES:  17,94.3 

Data  may  be  deleted  merely  by  entering  the  increment  of  that  data  point: 

ENTER  INTEGER  OF  DATA  TO  BE  DELETED:  7 

After  the  data  array  is  accepted,  you  will  have  the  option  of  printing  the  data.  Normally,  this  is  not 
required  as  the  array  will  be  printed  on  the  final  output. 

DO  YOU  WANT  TO  PRINT  THE  RAW  DATA?  N 

The  order  of  the  regression  must  be  entered.  Although  the  routine  can  process  ten  order  polyno¬ 
mials,  it  is  rarely  necessary  to  use  more  than  five  or  six. 

ENTER  ORDER  OF  LINEAR  REGRESSION  (LINEAR  =  1,  ETC):  4 

The  program  will  print  of  the  results  of  the  analysis  upon  completion  of  the  computations.  You  can 
print  the  analysis  if  desired: 

DO  YOU  WANT  A  COPY  OF  THESE  COMPUTATIONS?  Y 

Using  the  values  provided  in  figure  18,  the  resulting  polynomial  is: 

Y(i)  =  127.0306  -  9.31768*X(i)  -I-  1 . 178889*X(i)2  .  .  . 

where  Y:  longitude  and  X:  latitude. 

A  simple  x-y  plot  is  available  to  compare  the  fit  of  the  polynomial  to  the  raw  data.  If  selected,  you 
will  be  asked  to  enter  plot  title  and  axes  set  up.  You  will  have  an  opportunity  to  review  and,  if  neces¬ 
sary,  modify  the  inputs. 

DO  YOU  WANT  GRAPHIC  INPUT?  Y 

ENTER  HORIZ  (X)  AXIS  LEGEND:  LATITUDE 

ENTER  VERT  (Y)  AXIS  LEGEND:  LONGITUDE 

TITLE  PLOT  IS  E  ANDAMAN  SEA. 

HORIZONTAL  AXIS  IS  4  TO  16  WITH  TICK  MARKS  EVERY  1  UNITS. 
AXIS  LEGEND  IS  LATITUDE. 

VERTICAL  AXIS  IS  97  TO  103  WITH  TICK  MARKS  EVERY  1  UNITS. 
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AXIS  LEGEND  IS  LONGITUDE. 

IS  PLOT  SETUP  OK?  N 

SELECT  PARAMETER  TO  BE  CHANGED: 

(1)  PLOT  TITLE  (6)  MIN  VALUE  OF  VERT  AXIS 

(2)  MIN  VALUE  OF  HOR  AXIS  (7)  MAX  VALUE  OF  VERT  AXIS 

(3)  MAX  VALUE  OF  HOR  AXIS  (8)  TICK  INTRVL  OF  VERT  AXIS 

(4)  TICK  INTRVL  OF  HOR  AXIS  (9)  VERT  AXIS  LEGEND 

(5)  HOR  AXIS  LEGEND 

ENTER  NUMBER  OF  PARAMETER  TO  BE  CHANGED:  1 

The  analysis  may  be  repeated  without  entering  the  original  data: 

DO  YOU  WANT  TO  TRY  A  DIFFERENT  ORDER?  N 

You  now  have  the  option  of  recycling  the  program  with  different  data  or  ending  the  program. 

DO  YOU  HAVE  ANOTHER  DATA  SET?  N 

B2.  Routine  XYPLOT 

Start  Routine  XYPLOT  by  entering  XYPLOT.  You  will  be  asked  to  enter  the  plot  title: 

ENTER  PLOT  TITLE:  DEEP  SOCAL  BIGHT  —  PCNT  ZERO  SLD 

Several  plot  formats  are  available,  all  of  which  plot  statistical  data  as  a  function  of  month  using  print¬ 
outs  generated  by  Routine  XBT2A.  Options  1  through  3  plot  mean  and  standard  deviation  from  the 
mean,  while  the  fourth  option  plots  the  mean  only.  The  fifth  option  permits  you  to  plot  as  many  as 
six  different  parameters.  This  latter  option  is  helpful  in  illustrating  similarities  or  differences  among 
parameters.  All  plots  include  a  trend  line,  as  determined  from  the  January  and  December  means,  at 
the  start  and  end  of  the  monthly  data.  Typical  data  input  for  each  option  is  shown  in  brackets.  (Rou¬ 
tine  XYPLOT  does  not  include  this  information.) 

INDICATE  TYPE  UNITS: 

(1)  DEPTH  (MEAN  &  STD  DEV)  [layer  depth] 

(2)  TEMPERATURE  (MEAN  &  STD  DEV)  [water  temperature] 

(3)  FREQUENCY  (MEAN  &  STD  DEV)  [cut  off  frequency] 

(4)  PERCENT  (MEAN  ONLY)  [%  useful  SND  CHAN] 

(5)  OTHER  (USER  SELECTED)  [any  parameter] 

(6)  END  PROGRAM 
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ENTER  APPROPRIATE  NUMBER:  3 

The  assumption  is  made  that  all  input  data  are  in  metric  units.  However,  if  depth  or  temperature  data 
are  processed,  you  will  be  asked  to  select  output  units: 

INDICATE  UNITS  OF  OUTPUT  DATA: 

(1)  METRIC  (METERS,  DEGR  CELSIUS) 

(2)  ENGLISH  (FEET,  DEGR  FAHRENHEIT) 

ENTER  ( 1 )  OR  ( 2 ) :  1 

You  must  then  define  the  limits  of  the  vertical  axis  by  answering  the  appropriate  queries.  Note  that 
these  units  will  not  be  converted  to  English  units  if  English  units  were  selected.  In  most  cases,  the 
interval  between  tick  marks  will  be  determined  automatically. 

It  is  suggested  that  data  entry  be  made  using  the  keyboard  number  pad.  Note  that  the  keyboard  num¬ 
ber  pad  “plus”  sign  (+)  is  automatically  reconfigured  by  the  computer  as  a  comma  to  assist  in  data 
entry  for  Options  1  through  3  of  the  main  menu.  Leading  and  trailing  zeros  may  be  omitted. 
Depending  upon  the  option  selected,  the  computer  will  print  a  message  reminding  you  of  this  capa¬ 
bility. 

If  input  consists  of  mean  monthly  values  and  standard  deviation,  the  data  will  be  entered  in  the  fol¬ 
lowing  format: 

1  17.03+1.12 

2  16.54+0.1 

3  ... 

4  et  cetera 

If  standard  deviation  is  not  required,  only  the  mean  value  is  entered  at  the  prompt,  e.g.,  60.2.  In 
either  case,  you  are  given  the  opportunity  to  correct  errors  after  the  entire  array  has  been  entered. 

After  the  graphic  is  plotted,  you  will  have  several  seconds  to  examine  or  capture  the  screen  by  what¬ 
ever  method  you  select.  The  screen  will  be  erased  upon  pressing  <enter>.  You  will  then  be  given 
the  option  of  continuing  or  ending  the  program. 
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APPENDIX  C 

DEFINITIONS  OF  OCEANOGRAPHIC  TERMS 

Best  Depth:  The  best  depth  for  a  submarine  to  avoid  detection  by  a  hull-mounted  sonar.  This 
assumption  is  based  on  the  premise  that  maximum  downward  refraction  of  sonic  energy  will  occur  at 
the  depth  of  the  maximum  negative  sound  speed  gradient,  thus  causing  short  horizontal  sonar  ranges. 
For  the  purpose  of  this  paper,  the  depth  stratum  containing  the  minimum  gradient  must  have  a  thick¬ 
ness  of  10  m  or  more.  This  definition  disagrees  with  the  adage  that  defines  best  depth  to  avoid 
detection  as  sonic  layer  depth  plus  100  meters. 

Cut-off  Frequency  (COF):  The  maximum  frequency,  in  hertz,  that  will  remain  entrapped  within 
a  sound  channel.  Also  called  low-frequency  cutoff.  For  the  surface  layer,  COF  is  defined  as; 

Fo  =  0.3978  *  Voi-5  /  (SLD  *  (Vsid  -  Vo))®-^ 

where  Vq  and  Vsid  are  sound  speed  at  the  surface  and  sonic  layer  depth,  respectively.  Not  computed 
for  near  surface  layers  less  than  10  m  thick;  ignored  for  frequencies  greater  than  5.2  kHz. 

Cut-off  frequency  for  sound  channels  below  the  surface  duct  is  calculated  using: 

Fc  =  0.2652*Vai-5  /  (AZ  *  (Vb  -  Va))0-5 

where  AZ  is  sound  channel  thickness  and  Va  and  Vb  are  sound  speed  at  the  sound  channel  axis  and 
the  sound  channel  boundaries,  respectively.  Not  computed  for  sound  channels  less  than  10  m  thick. 

Deep  Sound  Channel:  A  subsurface  duct  for  long-range  transmission  of  acoustic  signals.  Axial 
depth  of  the  channel  is  the  depth  of  minimum  sound  speed  between  the  sea  surface  and  the  ocean 
floor.  Also  referred  to  as  the  primary  sound  channel. 

Half  Channel:  The  special  case  where  the  sound  speed  profile  exhibits  a  minimum  at  the  sea  sur¬ 
face  and  a  maximum  at  the  sea  floor.  By  definition,  neither  sonic  layer  nor  sound  channels  may 
occur  in  a  half  channel.  Normally  a  winter  phenomenon  associated  with  subpolar  regions  and  the 
Mediterranean  Sea. 

Near-surface  Layer:  The  layer  extending  from  the  sea  surface  to  sonic  layer  depth. 

Secondary  Sound  Channel:  A  sound  channel  located  between  sonic  layer  depth  and  the  axial 
depth  of  the  primary  sound  channel.  Also  referred  to  as  a  near-surface  sound  channel  to  distinguish 
it  from  the  deep  sound  channel.  Channel  thickness  is  the  depth  interval  between  the  sound  speed 
maxima  denoting  the  top  and  bottom  of  the  sound  channel.  More  than  one,  or  multiple  channels, 
may  occur  on  the  same  profile. 

Sigma-T:  An  abbreviated  form  of  density: 

Ot  =  1000  *  (q-1.0) 

where  q  is  density.  For  example,  sigma-t  of  28.5  equates  to  a  density  of  1.0285  gms/cms. 

Sonic  Layer  Depth  (SLD):  The  depth  of  maximum  sound  speed  above  the  deep  sound  channel 
axis. 

Surface  Channel:  An  acoustic  duct  between  the  sea  surface  and  the  sonic  layer  depth. 

Useful  Sound  Channel:  A  secondary  sound  channel  with  a  cut-off  frequency  of  less  than  800  Hz. 

Zero  layer  Depth:  The  special  case  where  sonic  layer  depth  occurs  at  the  surface.  Occurs  most 
frequently  in  tropical  and  subtropical  regions  where  surface  winds  are  insufficient  to  cause  mixing. 
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must  be  reformatted  into  regional  files  prior  to  analysis.  This  document  provides  guidelines  required  to  process  and  analyze  the  MOODS  data, 
thus  providing  valuable  information  required  to  characterize  the  oceanic  and  acoustic  properties  of  an  area. 
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